Predissociation of Diatomic Molecules: An Analytically Solvable Model 



Aniruddha Chakraborty 
School of Basic Sciences, Indian Institute of Technology Mandi, 
Mandi, Himachal Pradesh, 750001, India. 
(Dated: June 29, 2011) 

A model consisting of a Harmonic Oscillator well and a linear potential, coupled by Dirac delta 
function, is solved. We find the exact analytical expressions for Greens function for this problem. 
This Greens functions are used to calculate the Raman excitation profile and electronic absorption 
spectrum. 



INTRODUCTION 

O ; 

Nonadiabatic transition due to potential curve crossing is one of the most important mechanisms to effectively induce 
electronic transitions in collisions [H4l4j. Two state curve crossing can be classified into the following two cases 
according to the crossing scheme: (1) Landau-Zener (L.Z) case, in which the two diabatic potential curves have the 
same signs for the slopes and (2) non-adiabatic tunnelling (N.T) case, in which the diabatic curves have the o ppo site 
sign for slopes. There is also a non-crossing non-adiabatic transition called the Rosen-Zener-Demkov type [ll. |15|. in 
which two adiabatic potentials are in near resonance at large R. The theory of non-adiabatic transitions dates back 
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i— - n , to 1932, when the pioneering works for curve-crossing and non-crossing were published by Landau (16|, Zener [17 1 
• and Stueckelberg (l8j and by Rosen and Zener (l9| respectively. Since then numerous papers by many authors have 
been devoted to these subjects, especially to curve crossing problems[l], [l5| . In our earlier paper we have proposed 
an exactly solvable model for the two state curve crossing problem which assumes the coupling to be a Dirac delta 
function [9j. This model is used to calculate the effect of curve crossing on electronic absorption spectrum and on 
Resonance Raman excitation profile for the case of harmonic potentials [lfj. We have later generalized our model to 



deal with general multi-channel curve crossing problem too [111 ]. Even very recently our model ia extended to deal 
with nonadiabatic tunneling in an ideal one dimensional semi- infinite periodic potential systems 12j. Our work is in 
_^ , progress to deal with nonadiabatic tunneling in an ideal one dimensional finite periodic potential systems 13j . In the 
J> ■ present paper, we describe the predissociation of a diatomic molecule as a decaying quasi-bound vibrational state, to 
allow for an intersection of two electronic states. The model we propose to solve consists of two diabatic potential 
curves, approximated by a harmonic oscillator and a linear potential, coupled through a delta function at the crossing 
point. This model mimics the case in which we have an excited electronic state with a minimum, crossing a purely 



repulsive state. The advantage of this model is that it can be solved analytically. 

THE MODEL 

We consider two diabatic curves, crossing each other. There is a coupling between the two curves, which causes 
transitions from one curve to another. This transition would occur in the vicinity of the crossing point. In particular, 



it will occur in a narrow range of x, given by 

\V 1 (x)-V 2 (x)\^\V(x c )\. (1) 

where x denotes the nuclear coordinate and x c is the crossing point. V\ and V2 are the diabatic curves and V 
represent the coupling between them. In reality the transition between V± (x) and V2 (x) occur most effectively at the 
crossing, because the necessary energy transfer between the electronic and nuclear degrees of freedom is minimum 
there. Therefore it is interesting to analyze a model, where coupling is localized in space near x c rather than using a 
model where coupling is same everywhere (i.e. constant coupling). Thus we put 

V(x) = K 5(x-x c ), (2) 

here Kq is a constant. This model has the advantage that it can be exactly solved 0-[Hj]- 
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EXACT ANALYTICAL SOLUTION 



We start with a particle moving on any of the two diabatic curves. The problem is to calculate the probability that 
the particle will still be on any one of the diabatic curves after a time t. We write the probability amplitude as 

•<* 4 > -(*&))• (3) 

where if)\(x,t) and ?/> 2 (3;,i) are the probability amplitude for the two states. ^(x,t) obey the time dependent Schro 
dinger equation (we take h = 1 here and in subsequent calculations) 

dV(x,t) * , 

i V 1 > =H*(x,t). (4) 

H is dehned by 
where Hi(x) is 

H i (x) = -^—^ + V i (x). (6) 
Ira ox 

We find it convenient to define the half Fourier Transform ^(uS) of ^(t) by 



¥(«) = / (t)e""*di. (7) 
Jo 

Half Fourier transformation of Eq. (£[]) leads to 

= iG(w)*(0), (8) 

where G(cj) is defined by 

(w - iT)G(w) = 7. (9) 
In the position representation, the above equation may be written as 



where G(x,xq]lo) is 



Writing 



^{x,uj) = i I G(x,xo;uj) 1 $(xQ,ijj)dxo, (10) 



G(x,x ;uj) = (x\(w-H)- 1 \x ). (11) 



G(^ ;a,)=( ^f,-o;-) G 12 (,,, ;a;) , 
1 G 21 (a;,a;o;w) G 22 (a;, x ; w) ' 



and using the partitioning technique j20j we can write 

Gu(x,xq;uj) = (x\[u - H x -V(u> - H 2 )~ 1 V]- 1 \x ). (13) 

The above equation is true for any general V . This expression simplify considerably if V is a delta function located 
at x c 0411. 

n i s ra t \ i ^gG;(a;,a; e ;a;)G^(x c ,x c ;a;)G;(x e ,x ;a;) 

Gu( ^° ;W) = Gl(X ' X0 '" } + l-^G?(, Cl , c ;.)G0(, c ,x c ;c) ' (14) 
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FIG. 1: Schematic diabatic potential energy curves illustrating the model. 



where 

G° i (x,x ;u) = (x\(oj-H i )- 1 \x ), (15) 



and corresponds to propagation of the particle starting at xq on the second diabatic curve, in the absence of coupling 
to the first diabatic curve. Using the same procedure one can get [7H14j 

Similarly one can derive expressions for G 22 (x, xq;uj) and G 2 i(x, Xo; u>). Using these expressions for the Green's 
function in Eq. (j8]) we can calculate ^(uj) explicitly. 

The expressions that we have obtained for ^(w) are quite general and are valid for any Vi(x) and V 2 (x) (7rll3|. 
However, their utility is limited by the fact that one must know Gi(x,xq;lo) and G 2 (x, xq;u)). 



ELECTRONIC ABSORPTION SPECTRA AND RESONANCE RAMAN EXCITATION PROFILE : 

PREDISSOCIATION OF DIATOMIC MOLECULES 



In this section we apply our method to the problem involving Harmonic potential and linear potential. We consider a 
system of three potential energy curves, ground electronic state and two 'crossing' excited electronic states (electronic 
transition to one of them is assumed to be dipole forbidden and while it is allowed to the other) 0-Q • We calculate the 
electronic absorption spectra and resonance Raman excitation profile to understand the effect of predissociation. The 
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FIG. 2: Calculated electronic absorption spectra with coupling (solid line) and without coupling (dashed line). 



propagating wave functions on the excited state potential energy curves are given by solution of the time dependent 
Schrodinger equation 



,8_ ( ^f>(x,t) \ ( H mb ^{x) V 12 (x) 
dt\^f b {x,t) I V 21 (x) H vib>e2 {x) 



V>i (M) 



(17) 



In the above equation H v u,. e i(x) and H v n,. e2 (x) describes the nuclear motion of the system in the first electronic 
excited state (allowed) and second electronic excited state (forbidden) respectively 



and 



■6, el 



(x) 



lb.(2 



1 d 2 1 2 , , a 
2^dx^ + 2 mUJA(x - a) 



(x) 



1 d 2 
2m dx 2 



+ K F (x). 



(18) 



(19) 



In the above m is the oscillator's mass, lua is the vibrational frequency on the allowed states and Kp is the slope of 
linear potential for the forbidden states and x is the nuclear coordinate. Shifts of the nuclear coordinate minimum 
upon excitation is given by a, and V\ 2 (V 2 \) represent coupling between the Harmonic and linear potential, which is 
taken to be 

V 21 (x)=V 12 (x) = K 5{x-x c ), (20) 

where K represent the strength of the coupling. The intensity of electronic absorption spectra is given by [U, [22| 

I a (u)k Re[f? go dxf? ao dxo*r am (x) 
iG(x,x ;uJ + iT)W^ b (xo)}, 

where 

G(x,xq; uj + iT) = (x\[(uj q /2 + u - Lu eg ) + iT - ff u i6, e ] _1 |xo) . 

and 



Hvib.e — 



H vi b,el(x) K \X C )(X C \ 

Kq\x c )(x c \ H vibie2 (x) 



Here, T is a phenomenological damping constant which account for the life time effects. ^ tb (x, 0) is given by 







(21) 
(22) 

(23) 
(24) 
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FIG. 3: Calculated resonance Raman excitation profile for excitation from the ground vibrational state to the first excited 
vibrational state, with coupling (solid line) and without coupling (dashed line). 



where Xi( x ) is the ground vibrational state of the ground electronic state, luq is the vibrational frequency on the ground 
electronic state, ea is the energy difference between the excited (allowed) and ground electronic state. Similarly 
resonance Raman scattering intensity can be expressed in terms of Green's function and is given by 21, 22l |. 



iG{x, xq;uj - 



^dx ^f b *{x,Q) 

ir)*r b (zo,o)| 2 . 



(25) 



In the above i&f (x, 0) is given by 



y v ; b {x,o) 



x/O) 
o 



(26) 



where Xf( x ) is ^h.e final vibrational state of the ground electronic state. As G®(x, x ; oS) for the harmonic potential 
as well as linear potential is known 23J, we can calculate G(x, Xq; lu). 



Results using the model 

In the following we give results for the effect of curve crossing induced dissociation on electronic absorption spectrum 
and resonance Raman excitation profile in the case where one dipole allowed electronic state crosses with a dipole 
forbidden electronic state as in Fig. [TJ As in our earlier publications [7H9j, here the ground state potential energy 
curve is taken to be a harmonic one with its minimum at zero. The curve is constructed to be representative of the 
potential energy along a metal-ligand stretching coordinate [2l|. We take the mass as 35.4 amu and the vibrational 
wavenumber as 400 cm -1 Q for the ground state. The first diabatic excited state potential energy curve is displaced 
by 0.1 A and is taken to have a vibrational wavenumber of 400cto _1 . Transition to this state is allowed. The minimum 
of the potential energy curve is taken to be above 10700 cm -1 of that of the ground state curve. The second diabatic 
excited state potential energy curve is taken to be a linear potential. Transition to this state is assumed to be dipole 
forbidden. The two diabatic curves cross at an energy of 10804.1 cmT 1 with x c = 0.02477 A. Value of Kq we use in 
our calculation is Kq — 5.54275 x 10 -15 erg. A. The lifetime of both the excited states are taken to be 450 cm -1 . The 
calculated electronic absorption spectra is shown in Fig. [2] The profile shown by the dashed line is in the absence 
of any coupling to the second potential energy curve. The full line has the effect of coupling in it. The calculated 
resonance Raman excitation profile is shown in Fig. [3] The profile shown by the full line is calculated for the coupled 
potential energy curves. The profile shown by the dashed line is calculated for the uncoupled potential energy curves. 
It is seen that curve crossing effect can alter the absorption and Raman excitation profile significantly. However it is 
the Raman excitation profile that is more effected. 
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CONCLUSIONS 

In our earlier papers we have proposed an exactly solvable model for the two state curve crossing problem 
In this paper we have extended our model to deal with the case of predissociation. We have analyzed the effect of 
curve crossing on electronic absorption spectrum and on Resonance Raman excitation profile for the case of Harmonic 
potential coupled to linear potential. We find that the Raman excitation profile is affected much more by the crossing 
than the electronic absorption spectrum 0-0- As we have shown, our model allows for an analytic solution, with a 
simple physical interpretation. The diabatic representation is the natural one for our model. On the other hand, the 
almost constant coupling leads to an adiabatic representation. We would therefore be inclined to conclude that in 
those cases where the diabatic picture is more convenient, our model should be more reliable than those in which we 
only have a slight breakdown of the Born-Oppenheimer approximation. 
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